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Multiple Scattered Radiation Emerging from Continental Haze Layers. 

2: Ellipticity and Direction of Polarization 

By Gilbert N. Plass. George W. Kattawar, and Stephen J. Hitzfelder 

ABSTRACT 

The ellipticity and the direction of polarization are calculated for 
ladiation that has undergone multiple scattering from plane parallel layers. 
Both the radiation emerging from the top of the layer and that transmitted 
through the bottom are considered. Two different phase functions are 
used fcr the scattering layer: Rayleigh and haze L. The direction of 

polarization of the reflected radiation shows little variation as the 
optical depth of the layer increases, while there is a much larger varia- 
tion for the transmitted radiation. When the optical thickness is small, 
the direction of polarization for haze L varies rapidly with zenith angle 
near those angles at which the single scattered polarization is zero. The 
ellipticity of the radiation from haze L layers increases at first in 
direct proportion to the optical thickness of the layer. In general the 
ellipticity of the transmitted radiation is considerably greater than that 
of the reflected because of the greater average number of photon colli- 
sions in the former case. When the ellipticity is small, it is shown 
that the product of the polarization and the ellipticity is closely equal 
to -V/2I, where I and V are the first and fourth components of the Stokes 
vector . 


The authors are with the Physics Department, Texas A&M University, College 
Station, Texas 77843. 


I. 


Introduction 


The radiation scattered by a planetary atmosphere is polarized even 
though the incoming solar radiation is essentially unpolarized. In the 
first part of this article' (referred to hereafter as 1) we have discussed 
the degree of polarization together with the neutral points for photons 
emerging from the top and bottom of a homogeneous layer scattering either 
according to a Rayleigh or a haze L phase matrix. The dependence of the 
polarization on the optical thickness of the xayer and the solar zenith 
angle has been discussed in 1. 

Two other important parameters of the radiation field, the direction 
of the polarization and the ellipticity, are presented in this, the second 
part of this article. These results were ob'.lned at the same time and 
by the same matrix operator method discussed in I. They are presented 
separately for clarity and convenience of presentation. 

II. Direction of Polarization 

The angle X which the direction of polarization (maximum intensity 
component) makes with the direction of the meridian plane containing the 
final photon direction is given by 

tan 2 X - U/Q, (1) 

where U and Q are two of the four components of the Stokes vector (I,Q, 
U,V). The range 0° < X < 180° has been chosen for the principal value 
of X. The values of 0° and 180° have equivalent physical meanings. Tims 
when there is a continuous rotation of the d'r«**:t. 4 on -,-f polarization as 
u increases and X passes through 0° or L80°, this would be represented 
on a figure giving x as a function of P by a curve which goes off the 
top of the figure (180°) and reappears at the bottom (0°), or vice versa. 
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The direction of polarisation x of the photons scattered by a 
Rayleigh phase matrix from a layer of optical thickness t is shown in 
Fig. 1. Curves are given for both the photons reflected upward out of 
the top of the layer and for those transmitted through the layer and 
passing through its bottom surface in a downward direction. The curves 
in Fig. 1 are for the case when the cosine of the solar zenith angle 

b * 0.85332 (31.43°), the surface albedo A = 0, and the azimuthal angle 

o 

* 30° and 150°. The abscissa is the cosine of the nadir (for reflected 
photons) or zenith (for transmitted photons) angle of observation, V . The 
solar horizon is at the left of the figure, the antisolar horizon is at 
the right, and the nadir or zenith is at the center. 

For the reflected radiation and single scattered photons (i =0), 

X = 73° at the solar horizon and changes only a small amount as u in- 

O 

creases; near the nadir X decreases to 60 , a value that it must have as 
the nadir is approached from the 4> = 30° side for any optical thickness. 
The corresponding value of X is 120° as the nadir is approached from the 
i = 15C° side. As the angle of observation decreases from the nadir to 
the antisolar horizon, X varies rapidly at first and has the value 73° on 
the antisolar horizon. The range of X is taken as 0° to 180°; thus, if 
X is increasing on the graph and reaches the value 180°, the curve termi- 
nates at the top of the graph and reappears at the bottom at 0°. This 
is not a real discontinuity in the motion of the direction of polariza- 
tion, which moves in a perfectly continuous manner; only the graphical 
representation has this jump. A better representation is obtained if 
one imagines the top and bottom of the graph joined together to form a 
cylinder. Then there is no discontinuity in the curves at such a point. 
The variation of X with t is small for the reflected photons. The curve 
for T = 4 is a limiting curve whic!. is the same for all larger values of 
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t on the scale of the graph. 

The values of X shown in Fig. 1 for the transmitted photons exhibit 
more variation as the optical depth increases than was the case for the 
reflected photons. When the optical depth is very small, x decreases 
from a value of 73° at the solar horizon to 0° as <J> increases. Because 
of the chosen range of X, the value then jumps to 180° (as already ex- 
plained) and approaches 120° as the zenith is approached from this side. 
From the other side of the zenith (d> ■ 150°) X starts at 60°, increases 
to about 75°, and decreases to 73° at the anti-solar horizon. The curve 
for the photons transmitted through a layer of optical depth t » 4 shows 
variations of up to 20° from the curve from small t . However, the curve 
for t = 16 has an entirely different form. The direction of polariza- 
tion X has the value 90° everywhere except over a small region very near 
the zenith. The reason for this is that the photons transmitted through 
the diffusion region for layers of large optical depth have lost all 
memory of the original solar direction. 

The direction of polarization of radiation scattered by a Rayleigh 
phase matrix for a low solar elevation, U 0 * 0.18817, is shown in Fig. 2 
when <J> ■ 30° and 150°. The curve for t - 4 for the transmitted photons 
shows an interesting transition between the curves typical of smaller and 
larger values of t. 

The variation of x for radiation scattered by the haze L phase matrix 
is shown in Fig. 3. The variation is obviously more complicated than that 
for the Rayleigh phase matrix. There is little variation of y as the 
optical depth increases for the radiation reflected from the top of the 
haze L layer. On the other hand t depends strongly on the optical thick- 
ness of the layer for the radiation transmitted through a haze L layer 
of optical thickness t. When t “ 0 (single scattered photons), the varia- 
tion of x in the = 30° half plane is worthy of comment. The value of 


X is 150° at the zenith. As p decreases in the ■ 30° half plane, X 

decreases discontinuously by 90° as p passes through the value 0.9927 

2 

(9 = 6.93°). Kattawar and Plass have discussed the discontinuity 
which occurs whenever M (the second independent element of the scatter- 
ing matrix occurring in positions Sj and when the Stokes vector is 
chosen with the components I, Q, U, V) changes sign. Tor haze L the 
value of M changes sign at scattering angles (0) of 3.53° and 25.63°. 

When M - 0, the scattered radiation is unpolarized, and the value of X 
is indeterminate (see Fig. 18 of I for the zeros of the single scattered 
polarization). Thus there is no physical difficulty in determining the 
direction of polarization for any given scattering angle. 

When t * 0, P « 0.9927 and <J> * 30°, the scattering angle 0 3 25.63° 
and M “0; thus the value of x decreases by 90° as p passes through this 
value. As p decreases further in the d> * 30° half plane, the scattering 
angle decreases to 15.1° and then increases again passing through the 
value 25.63° again when l J ■ 0.6582 (6 ■ 48.83°). On^e more X decreases 
discontinuously by 90°. As p continues to decrease , X decreases con- 
tinuously until the horizon is reached. 

When t > 0, there is no actual discontinuity in the curve for x at 
the prints just discussed. For example, the curve .'or t ■ 1 follows fairly 
closely the t = 0 curve. There is a rapid chan*; u the value of X near 
P = 0.6582, but no actual discontinuity. When r 3 8, the curve is quite 
different from those for smaller optical thicknesses. When t ■ 16, x 
decreases slightly from the initial value of 150° at the zenith as p de- 
creases (<f = 30°) and then increases to a value of 180° at the horizon. 
Unfortunately the computer run was terminated at this value of the optical 
thickness. Presumably X would have the valve of 180° (or the equivalent 
0°) for all values o( W r - ue^* at the zenith for large optical thicknesses. 
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The variation of X in the <t> * 90° plane is shown in Fig. 4 when 

W “ 0.85332 (31.43°). There is little variation with optical depth for 

o 

the reflected photons, but considerably more variation for the transmitted 
photons . 

The variation of X for both the reflected and transmitted photons 
from haze L layers of various optical thicknesses is given in Fig. 5 for 

V * 0.188166 (79.15°). Once again there is little variation with optical 

o 

thicknesses for the reflected photons, but considerably more for the 
transmitted. There are no discontinuities in the X curves, since the 
scattering angle never becomes as small as 25.63° in this case. 

III. E lipticity 

The ellipticity of the radiation, E, is defined as the ratio b/a, 

where a and b are respectively proportional the major and minor axes 

of the ellipse described by the end point of the electric vector. This 

can be written in terms of the Stokes parameters 1, Q, U, V as 

_ 1 

E * - tan sin" 1 [V(Q 2 + IT + V 2 ) 2 ]} (2) 

and ^ 

E » - | V (Q 2 U 2 + V 2 ) 2 . (3) 

In almost all practical cases the ellipticity is small compared to unity 

and the approximate expression is sufficiently accurate. 

There is an interesting relationship between the polarization and the 

ellipticity. The polarization is defined as 

1 

P = (Q 2 + U 2 + V 2 ) 2 /I. (4) 

Thus 

E - - tan [~ sin -1 (V/PI) ] , 


(5) 


E . -V/2PI 


( 6 ) 


or 

2 I P E * - V (7) 

This simple approximate equation is usually quite accurate whenever E << 1. 

For haze L, Eq. (7) is accurate to six or seven significant figures for 
most angles. These equations remind us that the ellipticity is zero when- 
ever V is zero (which is always the case for Rayleigh scattering). The 
quantity V is also zero for any spherical polydispersion in the principal 
plane. 

The polarization in general tends to be large in a region where the 
ellipticity is small, as is shown by Eq. (7), and vice versa. It has been 
claimed that the calculated values of the ellipticity in the region where 
the polarization is small are uncertain due to truncation errors. We 
believe that our values are real from various tests which we have performed 
including varying the number of significant figures. All of our computed 
values of P, X (the angle of polarization), and E are very well behaved 
in the regions where E is small. 

The phase matrix for spherical particles has a zero entry in the 
fourth row and first column; thus incident unpolarized radiation cannot 
create a circular component from a single scattering. The photon must be 
scattered several times to be elliptically polarized. The ellipticity is 
shown in Fig. 6 for the reflected photons from a haze L layer when 

* 0.85332 (31.43°) and ■ 30° and 150°. A study of all the numerical 
results (including others not shown in the figure) shows that the ellip- 
ticity increases linearly with the optical thickness (t ) for small values 
of t . The ellipticity changes sign with V and thus is zero at certain 
particular angles. The absolute value of the ellipticity is plotted on a 
logarithmic scale in all the figures in this article. 
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The nadir angle at which the elllpticity is zero is given for certain 
solar angles in Table 1. In the <t> m 30° half plane, for example, there is 
one neutral point which does not change greatly with the optical thickness 
of the scattering layer for t < 1. When t - 2, there is no neutral point. 

On the other hand there are two neutral points when t - 4 (not shown in the 
figure). There is one neutral point for t > 8 near 0 • 16°. The be- 
havior of the neutral points for the elllpticity is Just as complicated 
as those for the polarization. 

The elllpticity of the reflected radiation does not increase appre- 

_3 

ciably with t for t > 0.25 and has typical values of about 10 away from 
neutral points. There may be zero, one, or two of these neutral points 
in each half plane in addition to the neutral point that always exists 
at the nadir. 

The ellipticity of the photons transmitted through a haze L layer 

of optical thickness t is given in Fig. 7 for the same solar zenith angle. 

Again the results (including others not shown in the fig re) show that 

the ellipticity increases linearly with t for small values of t. A few 

comments follow on some of the many complicated features shown in Fig. 7. 

A very sharp peak in the ellipticity develops near p = 0.66 in the $ = 30° 

half plane which persists for optical thicknesses as large as unity. There 

are neutral points whose position is given in Table 1. The ellipticity 

continues to increase in many regions as t increases up to t » 8 or 16. 

_2 

Values for the ellipticity between 10 and 0.13 develop in the $ ■ 30° half 
plane . 

It is important to note that typical values for the ellipticity of 
the transmitted radiation for larger optical thicknesses are from one to 
two orders of magnitude larger than the corresponding values for the re- 
flected radiation. The transmitted photons for the larger optical thick- 
nesses undergo more multiple scattering than the reflected photons. The 
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ellipticity increases in general with the number of scattering events, at 
least for relatively small numbers of such events. 

The ellipticity is always zero in the principal plane. The values when 
$ « 90° are shown In Fig. 8 as an example of the variation in another plane. 
The neutral points are given in Table 1. The other features already pointed 
out can also be seen here. 

The ellipticity when the sun is near the horizon (p * 0.188166, 

o 

9 ■ 79.15) is given in Figs. 9, 10, and 11 for the reflected and trans- 
mitted radiation in the <J> * 30° and 120° and the ■ 90° half planes. 

The linear dependence of the ellipticity on t when x is small is clearly 
seen in these figures. The neutral points are given in Table l. Their 
dependence on is especially interesting in the $ * 150° half plane for 
the transmitted photons. There is one neutral point when t < 0.0625; two 
when T ■ 0.25; none when * *= 1 and 2; and one when t >. 4. 

IV. Conclusions 

Ihe complicated variation of the direction of polarization and of 
the ellipticity of the radiatioi. scattered from the top and transmitted 
through the bottom of Rayleigh and haze L layers has been illustrated. The 
direction of polarization of the reflected radiation shows only relatively 
small variations as the optical depth of the layer is changed; on the 
other hand the variations are much larger for the transmitted radiation. 

The direction of polarization of the single scattered transmitted radia- 
tion for haze L may undergo a discontinuous change of 90° at those 
scattering angles where the polarization is zero. When the radiation is 
transmitted through a finite layer, the corresponding variation of the 
direction of polarization is rapid as the angle of view passes through the 
angle at which there is a discontinuity for single scattering (until the 
optical thickness becomes large). 
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The ellipticiCy of the radiation reflected by and transmitted through 
haze L layers Increases at first In direct proportion to the optical 
thickness (t) the layer as t increases from zero. In general the 
ellipticity of the transmitted radiation Is considerably greater than that 
of the reflected because of the greater average number of photon collisions 
in the former case. When the ellipticity Is small compared to unity, the 
relation V • - 2 1 P E if very accurate, where l and V are the first and 
fourth components respectively of the Stokes vector, P is the polariza- 
tion, and E is the ellipticity. 


This work was supported in part by Grant No. NCR 44-001-117 from the 
National Aeronautics and Space Administration. Acknowledgement is made to 
the National Center for Atmospheric Research, which is sponsored by the 
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Table 1 


Zenith r- nadir angle at which elliptlclty is zero 


Solar 

Zenith 

Angle 

Reflected (R) 
or 

Transmitted (T) 

Optical 

Thickness 

• 

o 

■ 

• 

o 

o 

• 

•e* 

4 • 150° 

31.43° 

R 

0.00195 

79.9° 

72.3° 

76.4° 



0.00781 

79.7° 

72.1° 

76.2° 



0.0625 

78.7* 

70.7° 

75.1" 



0.25 

76.9° 

67.9° 

72.4° 



1 

76.5° 

64.4° 

67.0° 



2 

None 

63.8° 

62.6° 



4 

6.3°; 13.3° 

64.0° 

58.7“ 



8 

16.3° 

64.0° 

58.0° 



16 

15.8° 

64.0° 

58.0° 

31.43° 

T 

0.00195 

75.5° 

72.2° 

78.2° 



0.00781 

77.3° 

71.9° 

O 

H 

• 

00 

f". 



0.0625 

76.2° 

70.4° 

77.1° 



0.25 

73.4° 

67.2° 

75.0° 



1 

64.9° 

62.4° 

72.8°; 86.6° 



2 

20.8° 

56.7°; 87.1° 

None 



4 

None 

69.9° 

44.1“ 



8 

None 

68.9° 

58.2° 



16 

None 

75.3° 

38.3° 

79.15° 

R 

0.00195 

None 

64.3° 

80.7° 



0.00781 

None 

63.0° 

80.7° 



0.0625 

None 

62.4° 

80.5° 



0.25 

None 

62.2° 

79.4° 



1 

None 

63.4° 

76.4° 



2 

None 

63.9° 

75.2° 



4 

None 

64.0° 

74.6° 



8 

None 

64.0° 

74.6° 



16 

None 

64.0° 

74.6° 

79.15° 

T 

0.00195 

88.9° 

64.1° 

77.5° 



0.00781 

88.5° 

64.0° 

77.3° 



0.0625 

85.7° 

63.1° 

75.9° 


Table 1 (cont.) 


Solar 

Zenith 

Angle 


Reflected (r) 
or 

Transmitted (T) 


Optical 

Thickness 

<p - 30° 

<J> » 90° 

4> - 150 

0.25 

80.2° 

61.0° 

71.2°; 85.4 

1 

63.1° 

54. 3-* 

None 

2 

38.3° 

42.2°; 79.5° 

None 

4 

None 

68.0° 

26.0° 

8 

None 

68.8° 

19.6° 

16 

None 

76.8° 

9.9° 


-14- 


Flg. 1. 

Fig. 2. 
Fig. 3. 
Fig. 4. 
Fig. 5. 
Fig. 6. 

Fig. 7. 

Fig. 8. 

Fig. 9. 


CAPTIONS FOR FIGURES 

Direction of polarization (X) as a function of cosine (p) of nadir 

angle (for reflected photons) and of zenith angle (for transmitted 

photons). The cosine of the solar zenith angle, P - 0.85332, the 

G 

surface albedo, A ■ 0, and the azimuthal angle, $ ■ 30° for the left 
half of the figure and 4> ■ 150° for the right half. Curves are given 
for various values of the optical thickness t of the Rayleigh scat- 
tering layer. 

Direction of polarization for u q * 0.18817, A * 0, $ • 30° and 150°, 

and Rayleigh scattering. See caption to Fig. 1. 

Direction of polarization for P q * 0.85332, A = 0, <f> * 30° and 150°, 

and haze L scattering. See caption to Fig. 1. 

Direction of polarization for P q ■ 0.85332, A » 0, * 90°, and 

haze L scattering. See caption to Fig. 1. 

Direction of polarization for P^ - 0.188166, A ■ 0, $ ■ 30° and 
150°, and haze L scattering. See caption to Fig. 1. 

Absolute value of ellipticity (E) as a function of cosine of nadir 
angle (P) for P p * 0.85332, A * 0, ■ 30° and 150° for photons re- 

flected from a haze L scattering layer of optical thickness t. 
Absolute value of ellipticity (E) as a function of cosine of zenith 

angle (p) for P^ = 0.85332, A = 0, ♦ *30° and 150° for photons 

transmitted through a haze L scattering layer of optical thickness T . 
Absolute value of ellipticity (E) as a function of cosine of nadir 
or zenith angle (P) for P q » 0.85332, A = 0, <(> * 90° for photons re- 
flected from and transmitted through a haze L scattering layer of 
optical thickness t. 

Absolute value of ellipticity (E) as a function of cosine of nadir 
angle (p ) for P q = 0.188166, A * 0, $ - 30° and 150° for photons 
reflected from a haze L scattering layer of optical thickness t. 
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Fig. 10. Absolute value of ellipticity (E) as a function of cosine of zenith 
angle (u) for ■ 0.188166, A ■ 0, $ - 30° and 150° for photons 
transmitted through a haze L scattering layer of optical thickness t. 

Fig. 11. Absolute value of ellipticity (E) as a function of cosine of nadir 
or zenith angle (u) for U p - 0.188166, A ■* 0°, $ - 90° for photons 
reflected from and transmitted through a haze I. scattering layer of 
optical thickness t. 
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